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ABSTRACT 

We analyze X-ray spectra of heavily obscured (7V# > 10 24 cm -2 ) active galax- 
ies obtained with Chandra, concentrating on the iron Ka fluorescence line. We 
measure very large equivalent widths in most cases, up to 5 keV in the most ex- 
treme example. The geometry of an obscuring torus of material near the active 
galactic nucleus (AGN) determines the Fe emission, which we model as a function 



of torus opening angle, viewing angle, and optical depth. The starburst/AGN 
O | composite galaxies in this sample require small opening angles. Starburst/AGN 

composite galaxies in general therefore present few direct lines of sight to their 
central engines. These composite galaxies are common, and their large covering 
fractions and heavy obscuration effectively hide their intrinsically bright X-ray 
continua. While few distant obscured AGNs have been identified, we propose to 
exploit their signature large Fe Ka equivalent widths to find more examples in 
X-ray surveys. 
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1. Iron Lines in the X-ray Spectra of AGN 

The strongest line in the X-ray spectrum of an active galactic nucleus (AGN) at mod- 
erate energies (4 < E < 10 keV) is due to iron fluorescence, particularly "neutral" Fe Ka 
at 6.4 keV, from Fe less ionized than Fe XVII. If both the continuum and the Fe-emitting 
region are viewed directly (the "Type 1" view), the equivalent width (EW) is small, typically 
less than 200 eV. As Krolik & Kallman (1987) pointed out, the EW can increase greatly 
if the fluorescing material is exposed to a stronger continuum than the observer detects. 
In this case (the "Type 2" view), an obscuring "torus" of material near the active nucleus 
blocks direct views of the central engine along the line of sight. X-ray observations of Seyfert 
galaxies generally support this unification scenario. Seyfert Is typically exhibit EW ~ 150 
eV (Nandra & Pounds 1994), while the EWs in Seyfert 2s reported hitherto tend to be larger 
and distributed more broadly, ranging from about 100 eV to 1 keV (Turner et al. 1997). 

Seen from the "Type 2" view, the Ka EW depends strongly on the torus geometry 
and total column density. Previous theoretical calculations (Awaki et al. 1991; Ghisellini, 
Haardt, & Matt 1994; Krolik, Madau, & Zycki 1994) have concentrated on producing EW 
< 1 keV, consistent with earlier observations. They indicated, however, that still larger EWs 
might result from Compton-thick tori (i.e., N H > 10 24 cm~ 2 toward the nucleus) with special 
geometries. 

Motivated by this suggestion and the recent discovery in Chandra X-ray Observatory 
(Chandra) observations of several much stronger lines (e.g., Sambruna et al. 2001), in this 
work, we examine the Fe Ka properties of active galaxies that have previously been identified 
as Compton thick. Our sample, listed in Table 1, comprises all such galaxies Chandra has 
observed for which published results, archival data, or proprietary data are available. Most 
of these are classified as Seyfert 2 galaxies. NGC 4945 lacks the requisite optical emission 
line signature of Seyfert galaxies, but X-rays reveal its active nucleus (Iwasawa et al. 1993). 
M51 is exceptional in this group for its low luminosity, and the source identified in the 
Chandra Deep Field South, CDF-S 202, is exceptional for its high luminosity, but they both 
fulfill the broad selection criteria. Although this sample is not complete and the data are 
heterogeneous, these measurements suggest that extremely large Fe equivalent widths are 
common in such heavily obscured AGN. 

2. Observations and Spectral Fitting 

All of the sample members were observed with the Advanced CCD Imaging Spectrometer 
(ACIS), and the Circinus galaxy and Mrk 3 were observed with the transmission gratings in 
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place. See Weisskopf et al. (2000), G. Garmire et al., in preparation, and Canizares et al. 
(2000) for more information on Chandra, ACIS, and the transmission gratings, respectively. 

Sako et al. (2000) analyze the Chandra spectra of Mrk 3. In order to examine the 
Fe lines in detail, we obtained these data from the Chandra Data Archive at the Chandra 
X-ray Center and extracted the High-Energy Transmission Grating Spectrometer spectrum 
following the method of Yaqoob et al. (2001), applying updated calibrations. Using XSPEC 
(Arnaud 1996), we fit the spectrum from 4-8 keV, avoiding the complication of many emission 
lines at softer energies while providing enough high-energy coverage to fit the H-like Fe and 
neutral Fe K/3 near 7 keV. We modeled the continuum as a strongly- absorbed power law, 
adding an unabsorbed power law constrained to have the same photon index to represent the 
scattered contribution. These data do not strongly constrain the intrinsic column density, 
but the EW is relatively insensitive to this uncertainty. (These data are in fact consistent 
with Compton thin obscuration; N H = 6.9(+2.4, —3.9) x 10 23 cm -2 .) Three Gaussians, which 
represent the Fe lines, are significant, with the neutral Ka EW= 0.71 keV. 

M51, NGC 4945, NGC 5135, and NGC 7130, were all observed with Chandra's back- 
illuminated S3 ACIS detector. We obtained the data on NGC 4945 and M51 from the 
Chandra archive. See Madejski et al., in preparation, for more complete analysis of the 
Chandra observation of NGC 4945. Terashima & Wilson (2001) present results from the 
observation of M 51, but we reanalyze the data here in order to utilize updated calibrations. 
The other two data sets are part of our program on starburst/AGN composite galaxies. We 
reprocessed all data from original Level 1 event files and applied the latest gain corrections, 
from 2001 July 31. We examined the lightcurves and excluded times during the NGC 4945 
observation where flares were significant. We detected no significant flares in the other 
data. In each case, we extracted the spectrum from a circular region of radius about 2" and 
measured the local background in nearby source-free regions. Except in the case of NGC 
4945, we did not bin the data, in order to retain the highest-energy photons. While the very 
soft X-ray spectra require thermal emission components, above 4 keV only the power-law 
continuum and Fe emission lines are significant. A single unabsorbed continuum component 
is required, which we interpret as the reprocessed contribution. The best-fitting continuum 
in each of these four spectra is relatively flat, consistent with pure Compton reflection. 

Table 1 summarizes these results, listing the measured central energy and the equivalent 
width of the neutral Fe Ka line for each sample member. We also record the continuum 
photon index, T, of the fits. The intrinsic continuum is strongly reprocessed in these obscured 
galaxies. While the simple power laws we apply are reasonable over the restricted energy 
ranges of the model fits, we caution that they are not appropriate descriptions over large 
energy ranges (e.g., for E = 0.5-30 keV). References in Table 1 include the identification of 
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the sample members' large column densities and published Chandra measurements, which 
we adopt for Circinus, NGC 1068, CDF-S 202, and IRAS 09104+4109. In these last two 
cases, Norman et al. (2002) and Iwasawa, Fabian, & Ettori (2001) apply reflection models, 
and the tabulated values of T represent the intrinsic source spectra. Thus, these two are 
steeper than the directly-measured continua of the other highly obscured examples, although 
their observed continua are similarly flat. 

3. Fe Line Models 

The Fe Ka EWs of the Compton thick AGNs in this sample tend to be extremely large, 
as expected when the intrinsic AGN continuum is fully hidden from direct view. We model 
the geometric dependence of continuum reprocessing with Monte Carlo simulations, which 
Krolik et al. (1994) describe more completely. Although we introduce no new physics in 
these calculations, the current observations demonstrate that more extreme configurations 
are relevant and must be evaluated. We consider a uniform torus of neutral material having 
a square cross-section of Thomson optical depth, r, in each direction, and opening angle 9 
(Figure 1). We consider all viewing angles, i, including Type 1 AGN orientations, where 
i < 9. We adopt intrinsic T = 1.9 and solar Fe abundance in the basic model and discuss 
additional variations below. We do not consider any scattered or direct contribution from 
ionized material that may be present outside the torus cavity or near the central source, but 
these simplifications do not significantly affect the EW. 

The EW rises smoothly with decreasing 9 and increasing i for constant r, as Figure 2 
illustrates for r = 4. For i < 9, the EW we calculate falls below 100 eV because we model 
only the torus contribution; the accretion disk, visible when i < 9, would add to the EW in 
that case. 

Figure 3 shows the variation of EW in the r — 9 plane for i = 65°. For fixed r, the EW 
is greatest when 9 is smallest, but EW does not increase monotonically with r. To achieve 
EW > 2 keV, r > 1 is required; the greatest EWs are found for r ~ 4. 

The case of 9 = 10° (Figure 4) demonstrates that extremely large EWs may arise over 
a wide range of viewing angles and optical depths for sufficiently small torus opening angle. 
While EW increases with i for a given r, moderate values of r and large values of i produce 
the largest EWs, which some of these Compton-thick AGNs require. 

These trends may be understood qualitatively as the interaction of several processes. 
Smaller opening angle means the torus obscures a larger solid angle around the nucleus, hence 
more Fe K-edge photons are captured and generate Ka photons. Viewing angles closer to the 
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equatorial plane see larger equivalent widths in the Ka line because the continuum is most 
severely attenuated along that direction. Large, but not too large, optical depth promotes 
Ka production by absorbing more hard X-rays: the continuum at Ka is suppressed, while 
photons well above the Fe K-edge are absorbed and trigger fluorescence. Compton optical 
depths more than ~ 4 diminish the equivalent width because the optical depth from the 
points of Ka production to the surface is too large for the fluorescence photons to escape. 

Variations in spectral shape, torus cross-section, and Fe abundance create changes in 
the EW that are only on the order of typical measurement errors. A shallower slope tends 
to increase the EW, because relatively more high-energy photons penetrate deeper into the 
torus, and the Fe-fluorescing region is then closer to the observer. Adopting T = 1.7, for 
example, the EW changes by less than 5% for % < 65°, and up to 30% at i — 75° for 9 = 30° 
and r = 2. 

Previous calculations have emphasized the effect of Fe abundance in order to account 
for large EWs, but abundance variations are only a secondary effect. With twice-solar 
abundances, for example, the EW increases up to 60% in the edge-on geometry. In contrast, 
purely geometric (and physically plausible) variations can alter the EW by an order of 
magnitude. 

The Fe line luminosity is correlated with the intrinsic AGN luminosity. For a given 
torus geometry, the Ka luminosity is relatively insensitive to viewing geometry, typically 
varying by factors of 5 over all viewing angles, with departures to factors of 10 only in a 
small minority of cases. 

We constrain the geometries of the obscuring regions of the galaxies with the largest 
EWs with these models. M51 requires 9 < 20°; for 9 = 10°, a range of optical depth r > 2 
and viewing angle i > 65° are consistent. NGC 5135 also favors 9 < 20°, allowing a range of 
% > 35°, depending on r. We note that ionized Fe lines are not detected in these two galaxies. 
Their large neutral covering fractions suggest that even if an ionized region is present, the 
neutral material may completely block it from our view. 

4. Significance 

Larger EWs require more extreme geometries, with torus opening angle 9 < 30° for 
EW> 2 keV. In this sample, the galaxies that exhibit the largest EWs have concentrated 
circumnuclear starbursts, while in the only galaxy certainly lacking a starburst, Mrk 3, the 
Ka EW is relatively small. This may be causal — the mechanical energy of a starburst, input 
as stellar winds and supernovae, may inflate the torus so that it covers a greater solid angle 
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around the nucleus (e.g., Wada & Norman 2002). 

The compact starbursts of NGC 5135 and NGC 7130 are evident from both their vacuum 
ultraviolet spectra, which show absorption features formed in the winds and photospheres 
of massive stars, and from their optical spectra where the high-order Balmer series and He 
I lines are observed in absorption (Gonzalez Delgado et al. 1998). The nuclear spectrum 
of M51 reveals high-order Balmer lines in absorption, also characteristic of a young stellar 
population (Heckman, Crane, & Balick 1980). The 100-pc scale starburst of NGC 4945 is 
detected in Paa images (Marconi et al. 2000), while a 200-pc ring of Ha and Br7 is evidence 
of the young stars around the nucleus of Circinus (Elmouttie et al. 1998; Maiolino et al. 
1998). The unpolarized nuclear spectrum of IRAS 09104+4109 exhibits a broad feature 
around He II A4686, attributable to Wolf-Rayet stars (Tran, Cohen, & Villar-Martin 2000). 
NGC 1068's starburst is large, and this extended (kpc-scale) star-formation region (Telesco 
& Decher 1988) therefore does not shape the immediate environment of the central engine. 
Instead, the AGN significantly affects the conditions of the nuclear region, as the prominent 
ionized Fe lines illustrate in this case (Ueno et al. 1994). 

Most recent investigations of the significance of starbursts in active galaxies have con- 
centrated on Type 2 AGNs (e.g., Gonzalez Delgado, Heckman, & Leitherer 2001) for the 
practical reason that the starburst signatures are more evident when the glare of the cen- 
tral engine is blocked. Our explanation for the exceptionally large Ka EWs seen in star- 
burst/AGN composites — very small opening angles — implies that any particular AGN with 
a starburst is more likely to be seen as Type 2, so few genuine starburst/Type 1 AGNs exist. 

Starburst/Type 2 AGN composite galaxies are common. Approximately half of all 
Seyfert 2s contain circumnuclear starbursts (Gonzalez Delgado et al. 2001). The composite 
galaxies are also preferentially more obscured than their "pure" counterparts, which lack 
starbursts (Levenson et al. 2001). The starbursts themselves are responsible for some of this 
obscuration. With typical star-formation rates per unit area SFR « 50-100M Q kpc~ 2 yr _1 
in the central 100 pc of composite galaxies (Gonzalez Delgado et al. 1998), the correlation of 
Kennicutt (1998) implies the mean gas column density in the star- forming regions is about 
10 24 cm~ 2 . The causal connection of large column density and star formation may then 
work both ways to enhance EW: strong starbursts arise where a large reservoir of material 
is available, and the starburst may also inflate the torus. Thus, in a significant fraction of 
AGNs, accretion occurs behind large obscuring column densities that cover most lines of 
sight. 

In order to reproduce the observed spectrum of the cosmic X-ray background, synthesis 
models (Setti & Woltjer 1989; Comastri et al. 1995) require a large population of obscured 
AGNs. An outstanding problem with these models has been the small number of obscured 
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luminous AGNs — Type 2 quasars — that are actually identified. Most of these are radio-loud 
sources, which represent only a minority of AGNs (Urry & Padovani 1995). Broad-band 
X-ray surveys are not generally an effective method for finding Type 2 quasars because even 
the X-ray continuum may not penetrate the enshrouding gas. For example, CDF-S 202 is 
the only high-luminosity obscured AGN reported in the 1-Ms exposure of the Chandra Deep 
Field South, and only one has been identified in the 185-ks observation of the Lynx field 
(Stern et al. 2002). 

This sample suggests that prominent Fe Ka (EW > 1 keV) is a common feature of very 
obscured AGNs, and the characteristic large EW can be exploited to find more of them. 
Specifically, we propose searching for Type 2 quasars in continuum-subtracted narrow-band 
X-ray images. At high redshift, in particular, where the host galaxy contribution is negligible, 
the Fe line contains nearly all the X-ray flux. Within a narrow energy range, AE 300 
eV, some objects might become significant detections that would be undetectable when the 
background or sensitivity constraint of a wider band is included. 

The Fe Ka line is a valuable probe of buried AGNs and the material that hides them. 
Isolating this line emission both spectrally and spatially from additional diluting sources, 
these Chandra observations suggest that large EWs are common in obscured AGNs. The 
most extreme EWs require small torus opening angles and arise in the starburst/AGN com- 
posite galaxies. 

We thank Tahir Yaqoob for essential assistance with and software for the grating spectra 
and thank David Strickland for providing a reduced spectrum of the NGC 4945 nucleus. 
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Table 1. Spectral Model Parameters 



Galaxy 




EW a 


r 


References 




(kcV) 


(kcV) 








M51 


6.43 ± 0.04 


4 9+ 215 


n 6 +2,6 

u -°-0.3 


1, 


2 


NGC 5135 


fi on+0.03 


24 +i.8 


+1 ' fi 
u.u_ 2 


3, 


2 


Circinus 


6.39 ±0.07 


2.11 ±0.48 


-0.8 


4, 


5 


NGC 4945 


6.39 + 0.02 


1 26+°- 26 


-0.98 ±0.7 


6, 


2 


NGC 1068 


6.40 


1.22 


l.otH 


7, 


8 


CDF-S 202 


6.43 + 0.34 


1 19+ 1 ' 2 


1.8 b 


9 


IRAS 09104+4109 


6.40 ±0.07 


1 1+0.9 


2.0 b 


10. 


11 


NGC 7130 


6 4n +0.05 


i i+0.3 
i - 1 -0.5 


-0 1+ ' 6 


12, 


2 


Mrk 3 


6.39 ±0.01 


n 7 o+0.46 
u -'°-0.37 


1.9 ±0.8 


13, 


2 



a Source rest frame 

b Reflection model fixed intrinsic source 

Note. — Errors are 90% confidence limits for one interesting parameter. 

References. — (1) Makishima et al. (1990); (2) this work; (3) Turner et al. 
(1997); (4) Matt et al. (1996); (5) Sambruna ct al. (2001); (6) Iwasawa et al. 
(1993); (7) Koyama ct al. (1989); (8) Young ct al. (2001); (9) Norman ct al. 
(2002); (10) Fraceschini et al. (2000); (11) Iwasawa et al. (2001); (12) Risaliti 
ct al. (1999); (13) Awaki et al. (1990) 
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Fig. 1. — Cartoon of model geometry. The square torus cross-section is shaded, and the solid 
line marks the torus symmetry axis. The torus opening angle, 9, and the viewing angle, i, 
are indicated along dashed and dot-dashed lines, respectively. 
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Fig. 2. — EW as a function of 9 and i, for r = 4 and T = 1.9. The image is scaled linearly 
from (white) to 10 keV (black). The dotted contour is at 100 eV, and the solid contours are 
scaled linearly in steps of 1 keV, beginning at 1 keV. Both 9 and i are indicated in degrees. 
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Fig. 4. — EW as a function of r and z, for 9 = 10°. The image and contour scales are the 
same as in Fig. 2. 



